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Separable Kernel of Nucleon-Nucleon Interaction in the Bethe-Salpeter Approach 
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The dispersion relations for nucleon-nucleon (NN) T-matrix in the framework of the Bethe- 
Salpeter equation for a two spin-one-half particle system and with a separable kernel of inter- 
action are considered. The developed expressions are applied for construction of the separable 
kernel of interaction for S partial-waves in singlet and triplet channels. We calculate the low en- 
ergy scattering parameters, the phase shifts, and the deuteron binding energy with the separable 
interaction. The approach can be easily extended to higher partial-waves for iVW-scattering 
and other reactions (NN-, nN- scattering). 
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Recent and planned experiments on TViV-scattering Q as well as reactions with the deuteron and light nuclei reach 
large momentum transfer from the probe particle to the investigated system. In this context, relativistic effects play an 
essential role in understanding reaction mechanisms. Studies of elastic electron-deuteron scattering (see, for example, 
|2|, |3[ ) in the nonrelativistic framework require to take into account mesonic exchange currents (a part of them has a 
pure relativistic origin, i.e. the so-called pair currents M). Theoretical calculations should therefore follow relativistic 
f~) ', approaches. 

The approach based on the Bethe-Salpeter equation (BS approach) J5| gives a powerful method to investigate 
various nuclear reactions both for electromagnetic and strong interactions. To find an analytic solution for the BS 
equation for T-matrix of NN scattering describing bound and scattering states, the separable kernel of interaction is 
very useful [pi . The separable kernel is well known from nonrelativistic consideration to solve the Lippman-Schwinger 
equation analytically M. The separable form of interaction is good not only for the technical reason but, in our 
.—, • opinion, also for the reflection of the nonlocal type of elementary nucleon-nucleon interaction (it was stressed first in 
papers H). This allows us not only to avoid the divergence peculiar of local theories but also to take into account 
f\ ■ the internal structure of the nucleon. 

In the paper, we consider the BS equation for T-matrix and solve it using the separable kernel of interaction for 

a spin-one- half particle system (for a spin-zero case see, for example, H and for trinucleon systems 0). This allows 

us to solve the BS equation without referring to the ladder approximation. We also find dispersion relations for the 

• i-H . obtained solution which allow us to perform explicit analytic calculations and connect internal parameters of the 

^^ ' kernel with experimental data: low-energy parameters, bound state energy, and phase shifts. 

The paper is organized as follows: after the partial- wave decomposition of the BS equation (section O), the one-rank 



separable kernel of interaction is introduced and used for solving the BS equation (section [II). Dispersion relations 
for the T-matrix are discussed in section IV . In section [V| it is shown that the separable kernel of interaction includes 



the kernel in the ladder approximation, and the connection between parameters of two kernels is found. Section VI 
is devoted to results of calculation of the low-energy parameters, bound state energy, and phase shifts of elastic 
TV TV-scattering and to the discussion for the future works. 
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II. PARTIAL- WAVE DECOMPOSITION OF THE BS EQUATION 

The BS equation for the two-nucleon T-matrix with momenta p\ (p[ ) , p 2 (p' 2 ) could be represented in the following 
form: 

T(p',p;P) = V( P ',p;P) + j^JdkV(p',k;P)S(k;P)T(k,p;P), (1) 

where P — pi + p 2 , p = (jp\ — P2)/2 [p' = (p[ — p' 2 )/2] are total and relative momenta, and V(p',p;P) is the NN 
kernel of interaction. The two-body propagator in terms of the total and relative momenta has the form: 

S(k; P) = 1/ UP/2 + k-m + ie) (P/2 -k-m + ie)l , (2) 

with a = a^ ■ 7'', and 7' 1 are Dirac matrices. 

After the partial wave decomposition in the rest frame of the two-nucleon system, the equation (|J) becomes (s = P 2 ), 

T a p(p' ,p',Po,p;s) = V a p(p' ,p',p ,p;s) (3) 

+ 7-3 5Z / dk ° / fc2 dkV ai (p' ,p' , k , k; s)S^g(k , k; s)T a p(k , k,p ,p; s). 

7(5 

(Here and below p = p\, p' = |p'|, k = \k\). Greek letters denote partial waves in the channels under consideration. 
In LSJp-basis |l(J, [H]) the singlet 1 S'o channel consists of four partial waves, and the triplet 3 Si channel, of eight 
partial-waves. 

Below we omit all states which have at least one negative-energy nucleon and take only ^(T-wave in the singlet 
channel and neglect also 3 D] t "-wave and take only 3 S^-wave in the triplet channel. Therefore eq. (^) could be rewritten 
as: 

t{p' Q ,p',Po,P;s) = v(p / ,p',p ,p;s) (4) 

i f„ f fc2& w (Po>P'> fc o,fc;s)£(fco,fc,.Po,.P;s) 



4tt 3 J u J (VS/2 -e k + ie) 2 - k 2 ' 

where t = Isa,?; = Vas, with a - 1 Sq or 3 S^ , and e k = Vk 2 + m 2 and m is the nucleon mass. 

The T-matrix normalization condition could be written in the on-mass-shell form. Introducing phase shifts <5(s), 
one could write 

t(s) ee t(0,p,0,P,s) = -— ^L==e lS ^ smd(s), (5) 

ysv s — 4m z 



with p = \/ s/A — m 2 = \/2mTi a b. 

To introduce also the low-energy parameters of iViV-scattering, it is suitable to expand the above expression into a 
series of p-terms with the following expression ||12f : 

pcotS(s) = -a^ 1 + r j-p 2 + 0(p 3 ). (6) 

Taking into consideration only the first two terms of the decomposition (ph, we define the scattering length a and 
the effective radius r$. 

As for the mass of the bound state (Mb), the squared T-matrix has a simple pole in the total momentum squared 
(s), and the bound state conditions could be written in the following form: 

j „' „-,„ „■ o _ ]\J2\ 



B(p' ,p / , Po ,p- ll =M 2 



t(P ,P ,Po,P,s) = —^2 1- R(Po,P ,Po,p;s), (7) 



where functions B and R are regular at the point s = M 2 . The bound state energy Eb is connected with Mb 
M b = 2m -E h . 



III. RANK I SEPARABLE KERNEL OF INTERACTION 

The ansatz for the separable kernel of interaction with rank I has the form: 

v(p' ,p',Po,p; s) = Xg(p' Q ,p')g(p ,p). (8) 

Here introduced is the g-function for NN- vertices, and A is the coupling parameter. 
To solve (|j) we could assume the following separable form for T-matrix: 

t(Po,p',Po,P; s) = T(s)g(p' ,p')g(p Q ,p) (9) 

and substituting (g) and (Q) into (||), we obtain the function t(s): 

r( S ) = l/(X- 1 + h(s)), (10) 

where the function h(s) is: 

^) = -i fdko fedk tbfo'W - (11) 



4^ 3 J U J (VS/2 - e fc + ie) 2 - fc 2 

Therefore, the T-matrix could be rewritten in the following form: 



t{Po,P,Po,p;s)= A _ 1 + ws i (12) 



and the on-mass-shell expression is: 



f(s) = d(i) = A-m,M ' ( } 

It should be noted that the expression (|l|) has the so-called N/D-iorm widely used in the nonrelativistic T-matrix 
theory p2] a nd in some relativistic methods. 

Using (03), it is easy to connect the solution of T-matrix and the phase shifts S(s). To achieve this, we assume that 
the imaginary part of the function n(s) satisfies the following condition: 

Imn(s) = 0. (14) 

This condition is connected with the specific choice of g-functions for NN- vertices to be discussed in the next section. 
By using ( [l3| ) and (|l4|), the phase shifts 5(s) could be expressed by the formula: 

co t S(s) = Rct(5) = - A ' 1 + Re Mg) (15) 

Imi(s) Im h(s) 

To express the low-energy parameters in terms of the T-matrix solution, it is suitable to expand the function h(s) 
in a series of p terms: 

h(s) = h + iphi +p 2 h 2 + ip 3 h 3 + 0{p 4 ), (16) 

Reh(s) = h +p 2 h 2 + O(p i ), (17) 

Im h(s)=p(h 1 +p 2 h 3 + 0(p 3 )). (18) 



Using now the definition (|6j) and ( |15[]l8| ), we obtain the ao and r$ parameters 

hi 2 



ao = A^TV r ° = 7^ 



(A- 1 +h Q )^-h 2 
hi 



(19) 



The bound state condition (0) with the help of ( |12| ) is now expressed in the form: 

A" 1 = -h{s = M 2 ). (20) 



IV. DISPERSION RELATIONS 

Let us consider now the analytic properties of the BS solution (namely, the function h) . The simplest choice of the 
function g(po,p) is the Yamaguchi type km: 

9(po,p) = {p 2 o-p 2 -l3 2 +ie)- 1 . (21) 

In this case, the function h could be rewritten as 

h(s,0) = - -r^dp2 / dp / p 2 dp — — = % — — — (22) 

4tt 3 M J J Ws/2 - e p + lef - p* p* - ej + te 

with dg2 = d/d@ 2 and ep — \/p 2 + (3 2 . We have introduced the second argument (3 to extract the explicit dependence 
of the function h on this parameter. 

Analyzing equation (E2J) , we could conclude that there are four poles in the complex plane poi namely: 

P { o ] W = ^ ~ e p + ze p[, 2) (,s) = - ^ + e p - ? e (23) 

(3)/ \ , • (4)/ \ 

With the change of the variable s, the poles p and p move, and there could be a situation when two poles 
"pinch" real pq axes. It means, the function h(s) has at this s-point the leap and imaginary part. First points at 
which this condition is satisfied (branch points) could be found from the following equations: 

PoH*) = -PoH") => s = 4m 2 , (24) 

P i o ) ( S ) = -Po 2 \s) => Sl =4(m + /3) 2 . 

Summarizing the situation, one could say that the function h(s) has two cuts starting at points sq and s%, respec- 
tively, and therefore could be written in a dispersion form: 

+ OG 

*»./^, (25) 

J s — s — it 

4m 2 

p(s', 0) = 0(s' - 4m 2 )p el (s\ (3) + 9{s' - 4(m + f3) 2 )p m (s', 0) 

with two spectral functions p e i-i n (el stands for elastic and in for inelastic) which are connected with the imaginary 
parts as 

p(s',0) = -lmh(s',0) = 1 ^(h-h*). (26) 

To find spectral functions, one should perform po-integration in (03) which gives the following result: 



1 



Ks,/3) - ~^dp J p dp s/4 _ v - sep + m2 _p 2 + . e 



1 1 



\fs — 2e p + it 2ep 



(27) 



Taking into account the following symbolic equation: 

1 V 



^finS(x - xq), (28) 



x — xq ± ie x — xq 
it is easy to find the spectral functions: 

s' > 4m 2 (29) 



p e l(s', (3) = v / 7v / s'-4m 2 /(7r 2 (s' - 4m 2 + 4/? 2 ) 2 ) (30) 

s' > 4(m + 0) 2 (31) 

p in (s', /3) = -(64/5 6 + 16,3V - 192/3 4 m 2 - 2O0 2 s' 2 
+ 192/? 2 m 4 - 320 2 s'm 2 + 16m V + 3s' 3 - 64m 6 - 12mV 2 )/ 
(2ttV(s' - 4m 2 + 4/3 2 )Vs' - Mm + 0) 2 ^fs' - 4(m - 0) 2 ). 



To perform integration (|22|), it is suitable to introduce new variables: 

s 
Am 2 



A/3 = — , 
to 




9 4ro 2 


. 1 


9 = 1 

s 


= 1 — 


2 TO + /3 

f = = 


l + A/3 



A 



w 



/3 



777,-/9 1 — A/3 

In case when the following conditions are valid 

to > (3 > =*> 

4(to + /3) 2 > s > Am 2 

the introduced parameters are real and positive: 

1 > p 2 > 0, 
v 2 > 1 > 0, 



TO 2 _ ^2 1 _ A/ ' 

a - 4(to + /3) 2 _ t-fl + A^) 2 
s-4(to-/3) 2 " t-(l-A/3) 2 ' 



1 > A/j > 0, 



(l + A /5 ) 2 >t>l, 



w z >0, 

cr 2 >0. 



(32) 
(33) 

(34) 



(35) 
(36) 



(37) 



The condition (pq) means that the second (inelastic) imaginary part does not contribute to the function lmh(s) when 
the phase shifts are calculated in the region 4(to + /3) 2 > s > Am 2 , and therefore: 



Im h(s,@) = Im/i e ;(s,/3) = 



p{ P 2 ~ i){\ + w 2 y 



if 4(to + /3) 2 > s > Am 2 . 



4TO 2 7r(/9 2 + W 2 ) 2 

Performing integration (Ea), we obtain for the real part of function h(s) (the imaginary part is given by (pa)): 



Reh(s,P) = hei(8,/3)+hin(s,0), 

hel up) = p(P 2 -l)(l + ^ 2 ) 2 ln 

4to 2 7t 2 (/7 2 + w 2 ) 2 



(38) 



(39) 



P+l 



hin(s) 



p-\ 

(p 2 -l)(l+77J 2 ) 2 ( W 2 -p 2 ) 
4to 2 7T 2 (/9 2 + W 2 ) 2 W 

(p 2 -l)(l + w 2 ) 
Am 2, n 2 {p 2 + w 2 ) ' 
{1 + n 2 ){\ + v 2 ) 2 

' 32TO 2 7T 2 (7J 4 + Cr 2 ) 



arctan 1/w 



In 



V 2 -l 



(l + a 2 )(l + 7; 2 ) 



2\2 



l6m 2 TT 2 (v 2 — l) 2 7j(u 2 — cr 2 ) 2 
x (v e + 7j 4 cr 2 - 4tj 4 + 4?j 2 cr 2 - v 2 - cr 2 ) arctan 1/v 

(1 + <J 2 )(1 + v 2 ) 3 

16m 2 ir 2 a(v 4 + a 2 )(v 2 — a 2 ) 2r K 2 {v 2 — l) 2 
x (-7j 6 + 2v 6 a 2 - 3tjV + 3ctV - 2a 4 + a 6 ) arctan 1/cr 
(l + a 2 )(l + 7; 2 ) 2 

l6m 2 TT 2 (v 2 — 1)(t; 2 — c 2 ) 

To find also expressions for low-energy parameters, we should return to (|18[)-([l9|) and expand the function h(s) in 
a series of p terms: 



ho(P) 



(w + (1 + w 2 ) arctan 1/w) 



(1 + w 2 ) 

Am 2 TT 2 w 3 

5 + 20tj 2 + 5v 8 + 14?j 4 + 20?/ /3tj 2 + 1 . /3v 2 + I 

= arctan 1/ 



(40) 



4m 2 7r 2( u 2 _ !)3( 3u 2 + X ) 

(tj 2 + l) 4 



Am 2 TT 3 v(v 2 — l) 3 



arctan 1/u — 



8777 2 7T 2 



■In 



U 2 +3 

(tj 2 + l) 2 
4m 2 7r 2 (?j 2 - l) 2 



6 

W) = ~ 4 4?r 2 5 ( 3W + ( 3 + ^ aiXtan V"0 ( 41 ) 

4m 4 7r 2 (w 2 - l) 5 (v 2 + 3)(3w 2 + l) 2 v 
+304i> 14 + 2704i; 10 + 29u 16 + 1212t> 4 + 2704u 6 + 29) 
/3w 2 + l , l3v 2 + 1 

x V^T3" arctan /1 

( V 4 + 10w 2 + l)(« 2 + l) 4 

- arctan l/i; 



4m 4 ir 2 v(v 2 - 1) 
1 -In 



8m 4 7r 



4^2 



llw 8 + 52v 6 + 66u 4 + 52w 2 + ll)(w 2 + l) 5 
8m 4 7r 2 (v 2 - l) 4 (3w 2 + l)(v 2 + 3) 



(1 + u; 2 ) 2 
ft x (/3 = I 1 + w ^ (42) 

4m 3 7TM; 4 

(1 + W 2 ) 2 (4 + 3 W 2 ) 

W) = 8^W ( 43 ) 

At this point we have anexplicit analytic expressions which connect the parameters of the separable kernel of 
interaction (A and (1 [\p\) with observables: phase shifts S(s) (eq. (|l5|)), low-energy parameters ao and r$ (eq. (|19|)), 
and bound state energy (eq. (p0|)). 

V. SEPARABLE AND ONE-MESON EXCHANGE KERNELS OF INTERACTION 

In this section, we show that the so-called realistic kernel of interaction (meson-nucleon) in the ladder approxima- 
tion is included in the separable kernel of the Yamaguchi-type g-functions. To achieve this, we introduce a simple 
approximation for the ladder kernel: 

V(po,p; ko, k) -> V(p ,p; k Q , fc) = — — . (44) 

Using the expressions of the kernel for scalar- meson (sc) exchange pi , |13| : 

g 2 1 1 
V sc (p ,p;k Q ,k) = -^^^-J. [( e p e k + m 2 )Qo{z fi ) - pkQ^z^)] , (45) 

where z M = (p 2 + k 2 — (po — ko) 2 + fi 2 )/(2pk), p is the exchange meson mass, and Qi(z) are the Legendre functions 
of second kind, Qq(z) — 1/2 In (z + l)/(z — 1), Qi(z) — zQo(z) — 1, it is clear that 

V sc {p o ,p;0, 0) = lim V sc (p ,p; k ,k) = Op5 M (po,p), (46) 

fco — >0,/c— *o 



with 

Dp 



9sc ! e p + m 



p Ann 2 2e p ' 

3m(P0,p)-1/(^-P 2 -M 2 )- (47) 

Expressions for U sc (0,0; fco, fc) could be obtained from eqs. (|4q-[47|) by the following substitutions: po — > fco and 
p — > k. To make connection between parameters more clear, we perform the p/m-decomposition in function a p up to 
0(p 2 /m 2 ) term: 

9 2 sc 1 
K c (0,0;0,0) = -^4ht (48) 

47T 7T Z (JL 



Using now expression (E3), we could write: 



2 

u sc (po,p; fco, fc) = -?r(-) 2 g»(Po,p)g»(k 0l k). (49) 



Comparing the above expressions with the separable form of the kernel introduced by m), we could write the 
following connection between parameters: 

P = H, \ sc = ~^f. (50) 

Equations (|50| ) are valid also for the vector-meson (vc) exchange kernel with the substitution of g sc — ► g V c and /i 
by the vector-meson mass and could be derived from the following expression: 

V vc (p ,p; fc , k) = -^£-L_J_ [-2(2e p e fe + m 2 )Q (z M )] . (51) 

The pseudoscalar-meson (ps) exchange kernel of interaction has the form: 

g 2 \ 1 
V ps (po,p;k ,k) = --p £ — -T-r [-(e p e fc -m 2 )Qo(^) + pfcQi(^)] , (52) 

and for p, fe — > we could write: 

V ps (p ,p;k ,k) ~ p 2 k 2 g f _ l (p () ,p)g tl (k Q ,k). (53) 

In this case, the expression for V ps goes to zero for p, k — ► 0, and it is impossible to find a connection between 
parameters similar to (|5C|). 

We have obtained connections between parameters for various cases of the scalar and vector exchange, which are 
given in Rcf. |1^] . Results are given in Table |[ 

VI. RESULTS AND DISCUSSION 

At this moment, we could find internal parameters of the kernel of interaction (A, /3) to reproduce experimental 
values for low-energy parameters a^ p , r^ p for the singlet channel (^Sq) and a^j. p and bound state (deuteron) energy 
Ef p for the triplet channel ( 3 5i). 

In the case of ^o-channel, we use equation ( |l9|) with clq = a e ^ p to find A: 

X- 1 = (a e x s P r l hx(P) - h (p). (54) 

Inserting the obtained expression into equation ( |20|) with vq = r^, we find: 

2 



exp 



[(osrr%(/3)-ft2G9)]- (55) 



hi((3) 

Solving the nonlinear equation (p5[), we could find the value of f3 and, then, using expression (p4), the value of A. 
In the case of 3 S'i-channel, we obtain A from the bound state condition ( p0[ ) with E\, = E d p : 

\- 1 = -h(s = (M e / p ) 2 ,P), (56) 

where Md = 2m — Ed- Inserting the obtained expression into (|19|) with ao = a e ^ p ', we find: 

exp _ hi(fi) 

° M -h (p)-h( s = Mipy [bt) 

Solving the nonlinear eq. (p7|), we could find the value of j3 and, then, using expression (|5q ) - the value of A. 

To solve equations @ and @, we use FORTRAN code with the DZEROX subroutine (as a part of CERNLIB 
package). As a result, we find the following values for A and /?: 



for 1 Sq channel : 

A = -0.29425404 GeV~ 2 , (i = 0.22412880 GeV, 

(58) 



for 3 S\ channel 



-2 



A = -0.79271213 GeV , (i = 0.27160579 GeV, 



TABLE I: Connection between parameters of two kernels. 



J p 




)i (GeV) 


g J /4ir 


A (GeV 2 ) 


0+ 


NNS 


0.983 


0.64 


-.06265954891 


0+ 


NNa' 


0.550 


7.07 


-.2166930823 


1- 


NNp 


0.769 


0.43 


-.02576448048 


1" 


NNlu 


0.7826 


10.6 


-.6577877886 



TABLE II: The binding energy and low-energy parameters for singlet and triplet channels. 



ISo 


a 0a (Fm) 


ros (Fm) 




Calculated 
Experiment 


-23.748 
-23.748 ± 0.010 


2.75 

2.75 ± 0.05 




35i 


am (Fm) 


rot (Fm) 


E d (MeV) 


Calculated 
Experiment 


5.424 

5.424 ± 0.004 


1.774 

1.759 ± 0.005 


2.224644 

2.224644 ± 0.000046 



The results on the binding energy of the deuteron and the low-energy parameters are shown in Table 0. The 
experimental data are taken from [14]. 

The phase shifts S s (s) and 8t(s) calculated with the above parameters are presented in Fig. pi The experimental 



data are taken from 15 . As it is seen from the figure, even with the simplest choice of the separable kernel of 
interaction - rank I with only two parameters A and /?, we find the low-energy parameters of elastic AW scattering 
a s and r s in the singlet channel and a t and bound state (deuteron) energy Ed in the triplet channel with required 
accuracy and reproduce phase shifts in a region till T; & = 100 — 150 MeV. 

To summarize, we recall once again the basic points of our lines of thought and the approximations made in the 
analysis. We start with the Bethe-Salpeter equation for T-matrix and perform the partial-wave decomposition. At 
this point, we omit all states which have at least one negative-energy nucleon and consider only 1 Sq -wave in the 
singlet channel and neglect also 3 D^ -wave and left only 3 5'^-wave in the triplet channel. In order to find the explicit 
analytic solution for T-matrix, we introduce separable ansatz for the kernel of interaction. As a result, we find the 
N / D-iovTii for on-shell T-matrix. Taking the simplest choice of g-function in a Yamaguchi form for AW-vertex, we 
analyze the analytic structure of the obtained solution for T-matrix and write the dispersion representation of the T- 
matrix. Performing the integration in the dispersion formula, we obtain explicit analytic expressions which connect the 
parameters of the separable kernel of interaction (A and (3 [A/3]) with observables: phase shifts, low-energy parameters 
and bound state energy. This connection is used to calculate intrinsic parameters of the kernel of interaction and then 
to find phase shifts. 

We find also that the kernel of interaction in the ladder approximation is included into the separable kernel. To 
obtain this result, we use a simple approximation for the expressions in the one-meson exchange model. 

The main result of the paper is the dispersion form of T-matrix for elastic AW scattering obtained for the separable 
kernel of interaction which allows us to perform analytic calculations and to connect explicitly parameters of the 
kernel and observables. The discussed method could be technically expanded to the multi-rank separable kernels and 
more complex forms of g-function for NN-vertex to achieve required accuracy for bound state energy, low-energy 
parameters and phase shifts in a wide region of cm laboratory energy. The approach can be easily extended to higher 
partial- waves for AW-scattering and other reactions (NN-, irN- scattering). 
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